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Abstract. Damped Lya absorbers are believed to be associated with 
galactic disks. We show that gravitational lensing can therefore affect 
the statistics of these systems. First, the magnification bias due to lens- 
ing raises faint QSOs above a given magnitude threshold and thereby 
enhances the probability for observing damped absorption systems. Sec- 
ond, the bending of light rays from the source effectively limits the 
minimum impact parameter of the line-of-sight relative to the center 
of the absorber, thus providing an upper cut-off to the observed neutral- 
hydrogen (HI) column density. The ray-bending also reduces a possible 
obscuration of the QSO by dust. The combination of the lensing effects 
yields a pronounced peak in the observed abundance of absorbers with 
high column densities 2 X 10 21 cm -2 ) and low redshifts (^ a bs <J 1) 
in the spectra of bright QSOs (B <J 18 mag) with redshifts £qso J> 2. 
The inferred value of the cosmological density parameter of neutral hy- 
drogen, i?Hl 5 increases with increasing redshift and luminosity of the 
sources even if the true HI density remains constant. This trend re- 
sembles the observed evolution of i?Hl(-^)- Damped Lya absorbers with 
column densities ;> 10 21 cm -2 and redshifts 0.5 <J z^bs <J 1 are re- 
liable flags for lensed QSOs with a close pair of images separated by 
~ 0.3" (v c /220 kms -1 ) 2 , where v c is the rotational velocity of the lens. 
Detection of these gravitational-lensing signatures with the Hubble Space 
Telescope can be used to constrain the depth of the absorber potential- 
wells and the cosmological constant. 
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Universe — Galaxies: Evolution — Quasars: Absorption Lines 



1. Introduction 

Most of the neutral hydrogen probed by QSO absorption spectra is traced by damped Lya 
absorption lines, corresponding to HI column densities of 10 20 ~ 22 cm -2 (see, e.g. Wolfe 
1988; Lanzetta, Wolfe, & Turnshek 1995). The high HI content and the cosmological 
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abundance of damped Lya systems suggest that they are associated with galactic disks 
(Wolfe 1988). Direct evidence for this identification is indeed provided by a wealth of 
observational facts. First, recent images of QSOs exhibiting damped Lya systems in 
their spectra show related galaxies (cf. Steidel et al. 1994, 1995). Second, abundant 
heavy-element absorption lines at low ionization stages are frequently associated with 
damped Lya absorption systems (Turnshek et al. 1989; Wolfe et al. 1993; Pettini et 
al. 1994). Third, the velocity field traced by these metal lines relative to the damped 
Lya line is consistent with typical galactic rotation velocities, if the damped component is 
associated with a galactic disk while the metal-line systems are embedded in a surrounding 
halo (Lanzetta & Bowen 1992; Lu et al. 1993; Turnshek & Bohlin 1993; Lu & Wolfe 
1994). Fourth, observations of redshifted 21-cm absorption and emission indicate disk-like 
structures that extend across galactic dimensions (Briggs et al. 1989; Wolfe et al. 1992). 
Fifth, Faraday-rotation observations of QSOs with damped Lya systems are consistent 
with the existence of micro-Gauss magnetic fields in the absorbers, of the same magnitude 
as detected in nearby disk galaxies (Welter, Perry, & Kronberg 1984; Wolfe, Lanzetta, & 
Oren 1992; Perry, Watson, & Kronberg 1993; however note the caveat in Kronberg, Perry, 
& Zukowski 1992). Based on these clues, we adopt the view that damped Lya absorption 
arises from neutral hydrogen in ordinary disk galaxies. 

The characteristic HI density profiles of nearby spiral galaxies can be used to get a 
rough estimate of the impact parameter of lines-of-sight to QSOs relative to the centers 
of damped Lya absorbers. Typical HI column densities of ;> 2 X 10 20 cm -2 imply impact 
parameters <J 15/i _1 kpc (Broeils & van Woerden 1994). This characteristic separation 
is indeed confirmed by direct imaging (Steidel et al. 1994, 1995). Geometrically, the 
absorbing disk can be intersected through a variety of impact parameters, with some lines 
of sight passing very close to the center of the absorber. If the damped systems resemble 
nearby disk galaxies, then the effects of gravitational lensing are strong and inevitable for 
lines-of-sight passing less than a few kpc from the absorber center, where the HI column 
density is high. The statistics of high column-density systems would then be marked by a 
pronounced lensing signature. 

Gravitational lensing gives rise to two distinct effects. First, it magnifies the flux 
of the background QSO population and thus increases the number of sources above a 
given magnitude limit (see, e.g., Narayan & Wallington 1993; Schneider, Ehlers, & Falco 
1992; Schneider 1992). The magnification bias tends to artificially increase the fraction 
of QSOs which show damped Lya systems. Second, the bending of light rays effectively 
limits the minimum impact parameter of the line-of-sight towards the QSO relative to the 
absorber center, and therefore provides an upper cut-off to the column-density distribution 
of damped Lya absorbers. The strength of these lensing effects depends on the depth of the 
gravitational potential-well of the absorbers, and on the underlying cosmology. Within the 
framework of Friedmann-Lemaitre cosmological models, a possible non-zero cosmological 
constant has a strong influence on the statistics of lensing events (Turner 1990; Kochanek 
1992). 

The magnification bias is most pronounced for high column-density absorbers in the 
spectra of bright QSOs. Systems specified by these two requirements are particularly easy 
to observe, given a sufficiently large sample of QSOs. The search for the lensing signature 
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requires space-based observations because lensing is primarily effective for absorbers at 
redshifts z a b s ~ 0.5-1 and sources at £qso J> 2. Unfortunately, existing surveys of damped 
Lya systems using ground-based and IUE observations (see, e.g., Lanzetta et al. 1995, 
and references therein) do not cover simultaneously the redshift ranges for absorbers and 
sources which are optimal for lensing. However, it should be straightforward to extend 
these surveys and to design an optimal search strategy that will probe the lensing signatures 
predicted in this paper through future observations with the Hubble Space Telescope. If 
gravitational lensing indeed occurs in damped Lya absorption systems, one would expect 
the formation of multiple QSO images especially in high column-density systems. It is, 
however, not surprising that multiple images were not resolved in such systems in the past, 
because the characteristic image splitting provided by spiral galaxies is <J 0.3". Damped 
Lya absorption at a fairly low column density, ~ 7 X 10 19 cm -2 , was observed in the 
spectrum of both images of QSO 0957 + 561 (Turnshek & Bohlin 1993), but in that case 
the large angular separation of ~ 6" is caused by an intervening galaxy cluster. 

The magnification bias also affects estimates for the value of the cosmological density 
parameter of neutral hydrogen in damped Lya systems, i?ni- It has recently been argued 
that the observed decline of between z ~ 3.5 and the present time (White, Kinney, 
& Becker 1993; Bahcall et al. 1993; Lanzetta et al. 1995) is a natural consequence of the 
conversion of gas into stars in disk galaxies (Kauffmann & Chariot 1994; Lanzetta et al. 
1995). Here, we examine the effect of lensing on the estimate of i?Hl(-^)- Although high 
column-density systems are rare their contribution to the estimate of 1S dominant, 
and so the lensing effect on their apparent abundance can change the apparent value of 
i?Hl considerably. In particular, we show that the inferred value of should generically 
decline with decreasing source redshift even if the true HI density in the absorbers remains 
constant. This apparent evolution results from the change in the efficiency of lensing with 
redshift. 

Thomas & Webster (1990) used a simplified model to examine the influence of the 
magnification bias on the observed statistics of metal absorption lines in QSO spectra. 
They concluded that the redshift distribution of the absorbers is at most weakly affected 
by lensing. Our study focuses on damped Lya absorption systems which are more likely 
to arise at small impact parameters where lensing is expected to be stronger. Moreover, 
we address additional statistical properties of the absorbers which are more sensitive to 
gravitational lensing. 

The outline of this paper is as follows. In §2 we develop the formalism required to 
investigate the various signatures of gravitational lensing. We adopt the simplest model 
for damped Lya absorbers, assuming no cosmological evolution in their HI density profiles, 
which we model according to observed profiles in nearby disk galaxies. This simple model 
provides a good fit to the observed column-density distribution of the damped Lya systems. 
Moreover, the strongest lensing signature is expected for absorbers at low redshift (^ a bs ~ 
0.5), where evolutionary effects due to mergers or accretion should be weak (Toth & 
Ostriker 1992; Zaritsky 1995). The numerical results for the various lensing signatures are 
presented in §3. Finally, §4 summarizes our main conclusions. 
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2. Formalism 



2.1. Parameterization of the Absorber Population 

In the following subsections, we calculate the probability that intervening disk galaxies pro- 
duce damped Lya absorption and simultaneously magnify the background QSOs whose 
spectra show the absorption lines. We assume that the absorber population can be pa- 
rameterized by the Schechter luminosity function (Schechter 1976), 

$(£) d£ = exp(-£) d£ , (2.1) 

where £ = (L/L*) is the luminosity in units of the luminosity scale Z**, and where, for 
spiral galaxies (Marzke et al. 1994) 

= 1.5 x 10 -2 h 3 Mpc -3 , s = -0.81 . (2.2) 

As usual, h is the Hubble constant Hq in units of 100 kms -1 Mpc -1 . 



QSOs 




Figure 1.- Schematic illustration of the lensing geometry. QSOs in the source plane are either multiply or 
singly imaged, depending on their distance from the optical axis. The two corresponding regions in the 
source plane are separated by the dotted curve, marked y = 1. Light rays are magnified by the lens on 
their way to the observer. The magnification factors, yUi and yU2 , of the two images of the central QSO 
are indicated for further reference. The inserted figure displays the face-on hydrogen column density as a 
function of radius in the HI disk. 
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We further assume that each light ray which passes an intervening absorption system 
face-on within a radius Rn(N') from its center shows damped Lya absorption with a 
column density equal to or greater than N' . Following common practice (e.g., Lanzetta et 
al. 1991; Wolfe et al. 1986), we take Rn(N') to be proportional to the optical radius R } 
and the optical radius to be proportional to some power t of the luminosity of the system, 

Rr(N') = h(N') R* t , (2.3) 

where R* is the optical radius of an L* galaxy. Observations suggest the values (e.g., 
Holmberg 1975; Peterson, Strom, & Strom 1979; Broeils & van Woerden 1994) 

R* w 11 h~ x kpc , tK, 0.4 . (2.4) 

Measurements of the surface density of neutral hydrogen in local disk galaxies (e.g., Broeils 
& van Woerden 1994) suggest that the surface density of neutral hydrogen is well fitted 
by an exponential function of radius. Figure 8b of Broeils & van Woerden (1994) shows 
average HI profiles such that 



MA-) = f-21o g ( 2xi g cm _ 2 ) (2.5) 



describes reasonably well the average disk radius within which the neutral-hydrogen column 
density exceeds N' . Note that N' is the face-on column density of the neutral-hydrogen 
disk. If the disk is inclined by an angle 6 relative to the line-of-sight, the observed column 
density is N = N' / cos 8. We neglect the finite scale-height of the hydrogen disk, which 
is typically smaller by 1-2 orders of magnitude than the disk radius. The geometrical 
cross section a\ jya (N) of a disk inside of which the observed column density exceeds N is 
obtained by averaging over all inclination angles, 



l 

2 



*L ya (N) = *R' / d7/H(^7)7H( 7 -7o), (2-6) 
Jo 

where 7 = cos 6, and the Heaviside step function H(7 — 70) accounts for the fact that very 
high column densities can only be observed in disks which are inclined by an angle 6 > 80 , 
with 

VlO 3 / 4 



70 = cos #0 = niin 



2 X 10 20 cm- 2 



(2.7) 



N 

Evaluating the average in equation (2.6), we obtain 

^Lya(iV) = 7TR 2 (/|(iV)) , (2.8) 

where 



7o 



(fl(N)) 1/2 2 



1/2 

&(N) - 2/ H (iV) 21i ; n 70 io ~ 1 + ^(1 - 21n 7 o + 21n 2 70) 



(2.9) 



We define the angle- averaged disk radius for absorption with column density > N by 

(R R (N))=R(fl(N)) 1 / 2 . (2.10) 
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We model the lensing properties of the intervening systems as singular isothermal 
spheres, assuming that these systems resemble galactic disks with flat rotation curves 
which are embedded in spherical halos. Such systems are the simplest to consider, as they 
are characterized by only one parameter, viz. their one-dimensional velocity dispersion, 
o v . HST observations of bulges of spiral galaxies show no indication for cores, i.e., central 
regions in which the surface-brightness is constant (Tremaine et al. 1995, in preparation). 
Similarly, models for observed gravitational-lens systems indicate that the density profiles 
of the lenses are either singular or have very small core radii <J 200 pc (Rix, Schneider, & 
Bahcall 1992; Kochanek 1991). Assuming a circular rotational velocity of u C; * = 220 km s _1 
for an L* galaxy (e.g., Peebles 1993; Aaronson, Huchra, & Mould 1979), we get 

^ * = ^ ^ 160 kms" 1 . (2.11) 

The rotational velocity, and hence the velocity dispersion, are related to the luminosity 
through the Tully-Fisher relation (Tully & Fisher 1977), 

a v = I 1 '" a Vj * . (2.12) 

The Tully-Fisher exponent a depends on the waveband where the observations are taken; 
it ranges from a 2 in the blue to a 4 in the infrared (see Strauss & Willick 1995, and 
references therein). 



2.2. Probability for Damped Lya Absorption 

Ignoring gravitational lensing, the probability for a QSO to exhibit damped Lya absorp- 
tion in its spectrum can be calculated by dividing the physical path length surveyed for 
absorption by the mean free path between subsequent absorptions, 



^Lya(iV) = n Lya( r Lya (iV) AX 



(2.13) 



where nL ya is the spatial density of absorbers and a\ jya (N) is their individual cross sec- 
tion for damped Lya absorption with a column density > iV, given by equation (2.8). 
Integrating over galaxy luminosities we find 



n Ly a^Lya(iV) = 0* T(l + s + 2t) ttRI {&(N)) 



(2.14) 



where the Gamma function results from the Schechter luminosity function, taking the 
luminosity dependence of R into account. The physical path length between redshifts z\ 
and Z2 is given by the absorption-distance difference, 



AX(z u z 2 ) 



Z2 



dz (1 + zf 



■dt 



dz 



(2.15) 



where the factor (1 -\-z) 3 accounts for the change in the galaxy density due to the cosmolog- 
ical expansion. The proper-distance interval \cdt \ depends on the cosmological parameters, 



■dt\ 



c 



dz 



1 



i + * V n o(i + zf + (i-n - n A )(i + zf + n y 



(2.16) 
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where i?o is the present-day density parameter, and Qa is the cosmological constant. 
Combining the last four equations, we obtain 



P Lya ( Zl ,z 2 ,N) = $,r(l + s + 2t) 7rRl(fl(N)) I * dz (1 + zf 



■dt 



dz 



(2.17) 



(see also Bahcall & Peebles 1969, and Burbidge et al. 1977). Note that P\ jya (zi } z 2} N) is 
independent of the Hubble constant. 

This equation changes considerably when gravitational lensing is taken into account 
due to the magnification bias. The derivation of the absorption probability can then be 
split into two steps; one to calculate the probability for magnification, and the second to 
convolve the magnification probability with the QSO number counts. 

We adopt the conventional description of gravitational lensing (e.g., Schneider et al. 
1992). The optical axis of the lens system is defined by the positions of the observer and the 
center of the lens, and the lens- and source planes are constructed orthogonal to the optical 
axis at the redshifts z and z s of the lens and the source, respectively. Later, we shall refer 
to z and z s by the redshifts of the absorbers, £ a bs, and the QSOs, £qso- The lens plane 
constitutes the observer's sky in the vicinity of the lens. We normalize the coordinates in 
the lens plane by the length scale appropriate for singular isothermal spheres, 

£„ S 4,(^) J f ™ S 4»(^) J (2.18) 

V c / iio r s V c / iio 

where ra,r s , and ra s are the angular-diameter distances from the observer to the lens 
(deflector) and the source planes, and from the lens to the source plane, respectively, in 
units of the Hubble length (c/Ho). For a detailed description of the lensing properties of 
singular isothermal spheres, see §8.1 of Schneider et al. (1992). Coordinates in the source 
plane are correspondingly scaled by the projection of £o from the lens- onto the source 
plane. Using equation (2.12), the length scale £o can be written as 

e =4vr(^) 2 -Z-d(z,z a )P' a = -Z-0.d(z,z B )(?' a . (2.19) 

We denote the scaled coordinates in the lens and source planes by x and y, respectively. 
Then, sources with y < 1 have two images at x = y ± 1, while sources outside y = 1 have 
a single image at x = y + 1. The total magnification of a point source is 

a = i 2 ' y {y ^ l \ . (2.20) 

[ 1 + 1/y (otherwise) 

If a source is multiply imaged, and if the multiple images are not resolved when the 
spectrum of the source is taken, then the damped Lya absorption of the two superposed 
spectra is dominated by that of the brighter image which is formed at the larger impact 
parameter, x = y + 1. In order to get damped Lya absorption with column density > N 
in the composite spectrum, the outer image of the source must lie at a sufficiently small 
impact parameter. This is equivalent to the requirement 

y < m(N) = x K (N) - l = (i?H c (iV)) - 1 . (2.21) 

?0 
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Inserting equations (2.3), (2.10), and (2.19) into (2.21), we obtain 

VEi^N) = H '' R f} (N) r I'- 21 " ~ 1 • (2.22) 

cfl* d{z, z s ) 

For reasonable choices of t and a, (t — 2/a) ~ 0, so that yji(z } z S} N) depends only weakly 
on £. In the following, we shall assume that yn(z } z S} N) is independent of £ since this 
simplifies our discussion considerably. 

The dimensionless cross section of the intervening system for magnifying a point source 
by more than a factor /i and causing damped Lya absorption with a column density of 
> N is defined to be the corresponding area in the source plane within which a point 
source has to lie, normalized by the square of the length scale in the source plane. For the 
singular isothermal sphere, the cross section is given by (cf. Schneider et al. 1992) 

a(z,z s ,N;ri = { 47r lf ^ 
v ' ' \Try^(z,z s ,N) (otherwise) ' v ; 

if y K (z,z s ,N) < 1, or by 

( 47r//i 2 (/i > 2) 

<t(z,z s ,N;(j.) = I vr/(/i-l) 2 [2 > fi > 1 + l/ m (z,z s ,N)] (2.24) 
y 7ry|[(z, z s , N) (otherwise) 

if yn(z,z s ,N) > 1. 

The probability for a QSO to be magnified by more than /i and at the same time 
to show damped Lya absorption with column density > N in its spectrum can now be 
found by calculating the fraction of the total area of the source sphere covered by the 
appropriate cross sections of the lens population. Implicitly, this assumes that lenses and 
QSOs are physically uncorrelated and that the cross sections do not significantly overlap; 
both of these assumptions are well justified. Integrating over source luminosity and redshift 
between redshifts z\ and Z2 < z s , we obtain the probability that a QSO at redshift z s shows 
damped Lya absorption with column density > iV, and is magnified by a factor > /i by 
lenses in the redshift interval z\ < z < ^2, 

p' Glj (z 1 ,z 2 ,z s ,N ] tj)= ^elr(i + s + - 

W V a ' (2.25) 

f Z2 cdt\ 
X / dz d 2 (z,z s )(l + z f 



dz 



a(z } z S} N;/i) 



As in equation (2.14), the Gamma function results from integrating over luminosity, as- 
suming that the dimensionless magnification cross section is independent of £. 

Equation (2.17) should be a limiting case of equation (2.25) in the sense that if we 
disregard lensing, (2.25) should equal (2.17). In fact, if we allow for any magnification 
factor by setting /i — > 1, the cross section is limited by the column-density requirement 
only, and we obtain from equations (2.23) or (2.24) 

a(z,z s ,N;n) -> iry^z, z S} N) . (2.26) 
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If /i — > 1, we further have yji(z } z S} N) ^> 1, and thus 



<r(z,z s ,N;ri^n(&(N))Rl 



c6* d(z } z s ) 



(2.27) 



Inserting this into equation (2.25), we obtain 



P^ L (z 1 ,z 2 ,z s ,N; f i)^^r[l + s + -) rrRl (f*(N)) 



Z2 



dz (1 + zf 



■dt 



dz 



(2.28) 



which is identical with equation (2.17) provided that t ~ 2/a, as we have assumed before for 
simplicity. Note, however, that lensing reduces the area in the source plane within which 
QSOs show damped Lya absorption with column density > N because of the focusing 
effect of the lens. More precisely, this area is reduced compared to the case when lensing 
is ignored by the average magnification factor within yn(z } z S} N). 

Let now (dNQso / dS)(S) dS be the intrinsic number of QSOs at redshift z s with flux 
between S and S + dS. Due to the magnification bias, the observed number of QSOs with 
flux greater than S which show damped Lya absorption with column density > N between 
redshifts z\ and z 2 is given by 



JVqso(S) 



f°° ( S 

J dS' P l G l J [zl,z 2l z Sl N■ 1 — 



S \ dN "qsq 
dS 



(S>) 



(2.29) 



(see, e.g. Schneider et al. 1992; Narayan & Wallington 1993). Therefore, the probability 
for a QSO at redshift z s to be detected with a flux greater than S and to show damped 
Lya absorption with column density greater than N is given by 



If 00 / S 

Pgl(zi 7 z 2} z S} N; S) = — — / dS 1 P GL [ z 1} z 2} z S} N; — 

JVqso(^>) J \ b 



S\dN QSO{s , h (230) 



dS 



If we disregard gravitational lensing, Pgl(^i , z 2 , z s , N; S) reduces to equation (2.17) and 
becomes independent of the QSO flux S. 

The intrinsic QSO number counts are well modeled by a broken power law (e.g., Boyle, 
Shanks, & Peterson 1988, Hartwick & Schade 1990) 



dN< 



QSO 



dS 



(S) 



S~ a (S<S ) 
S~ b (otherwise) 



(2.31) 



where a = 1.64 and b = 3.52 provide a reasonable fit to observational data (Pei 1995). The 
break luminosity So corresponds to an apparent QSO magnitude of B 19.5. 
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2.3. Inferred Properties of Damped Lya Absorbers 

One of the goals of observing Lya absorbers is to infer the amount of neutral hydrogen 
in the universe, and thereby to study the history of systems which contain neutral hydro- 
gen. An intermediate step in this calculation is to derive the column-density distribution 
function f(N) of the absorbers, which is defined to be the number of absorbers with col- 
umn density within dN of N per unit absorption distance dX . Using the definition of the 
absorption distance, the probability without lensing (2.17) can be written as 



P Lya (^ , z 2 , N) = 0* T(l + s + 2t) vr Rl (&(N)) [X(z 2 ) - X( Zl )] 



(2.32) 



If among a sample of q QSOs / damped Lya absorbers are found within a column-density 
range of AN around iV, and if the absorption distance surveyed per QSO is AX(z\, z 2 ) } 
then the column-density distribution inferred from these observations is given by 



f(N) dNdX = — 



I 



H qAX{z u z 2 )AN 



dNdX 



(2.33) 



Since {I I q) approximates the probability for a QSO to show the specified absorption, we 
can write 



f(N) 



1 



dP{z u z 2 ,N) 



dN 



(2.34) 



Ho AX(z u z 2 ) 

Note that f(N) scales linearly with the value of the Hubble constant, because P is inde- 
pendent of Ho and AX is proportional to Hq 1 . Ignoring gravitational lensing, equation 
(2.32) can be substituted into equation (2.34) to give 



hya(N) 



H< 



■$*r(l + s + 2t)TrRl 



dN 



(2.35) 



This result is independent of the QSO redshift, the QSO flux, and the redshift interval 
surveyed for damped Lya absorption. If, however, gravitational lensing is taken into 
account, equation (2.34) reads 



fGh(zi,z 2 ,z s ,S;N) 



1 



Ho AX(z u z 2 ) 



dP G h(zi,z 2 ,z s ,N; S) 



dN 



(2.36) 



The resulting column-density distribution depends on the redshifts selected, the cosmo- 
logical parameters, and, in particular, on the QSO flux S because of the magnification 
bias. 

Given f(N) } the inferred comoving density of neutral hydrogen in systems with column 
densities of Ni < N < N 2} normalized by the present critical density, is found to be 



nm(N 1 ,N 2 ) = H ± m 



N 2 



dNNf(N) 



(2.37) 



C Pc,0 JN 1 

where rh is the mean molecular mass per proton, and p c q is the present-day closure density, 



Pc,0 



3HI 

8ttG 



(2.38) 
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The lower integration bound is fixed by the requirement that the absorption be damped, 
i.e. N\ ~ 10 20 cm -2 . In principle, the choice of N2 is arbitrary, but in existing surveys it 
is limited to N2 <J 10 22 cm -2 . Since f(N) scales linearly with the Hubble constant, 
is independent of Ho . 

The fact that f(N) depends on the redshifts chosen and the QSO flux when gravita- 
tional lensing is taken into account implies that the inferred value of depends on these 
parameters. 

3. Results 

Based on the formalism developed in §2, we now consider quantitatively the characteris- 
tic signatures imprinted by gravitational lensing upon quantities inferred from studies of 
damped Lya absorption systems. We begin with the column-density distribution func- 
tion, f(N). Figure 2 shows N f(N) for two choices for the absorber-redshift interval and 
for three cosmological models. We fix the QSO redshift at £qso = 3 and consider the 
absorber-redshift intervals of 2 < z a b s < 3 (Fig. 2a) and 0.4 < z a b s < 1 (Fig. 2b). The 
heavy line in both panels of figure 2 shows N f(N) ignoring gravitational lensing, and the 
three additional curves display the effect of lensing for the cosmological models i?o = 1, 
£2a = (solid line); i?o = 0.2, Qa = 0, (dotted line); and i?o = 0.2, Qa = 0.8 (dashed 
line). 

The observed column-density distribution of Lya absorbers can be well fitted by a 
power law over a broad range of column densities, f(N) oc N@ with fi = —1.7 ±0.2 (Tytler 
1987; Lanzetta et al. 1991). Our results for f(N) reproduce this power-law behavior well 
for N <J 10 21 cm -2 where the effect of lensing is weak. This indicates that our modeling 
of the neutral-hydrogen profiles in galactic disks is reasonable. 

For high-redshift absorbers (^ a bs <J zqso) 5 N f(N) drops below the result expected 
without lensing at N ~ 10 21 cm -2 , then rises to a peak around N ~ 10 22 cm -2 , and 
finally declines sharply at still larger N. As mentioned in §2, the area in the source 
plane within which QSOs show damped Lya absorption with a given column density is 
reduced compared to the geometric cross section by the focusing effect of gravitational 
lensing. The magnification bias compensates this reduction for increasing iV, because 
lines-of-sight passing through higher column densities experience higher magnifications on 
average. The magnification bias is weak for high-redshift absorbers (cf. Fig. 2a). However, 
for low-redshift absorbers and high-redshift QSOs, lensing by the absorbers is much more 
efficient, so that the peak in N f(N) } produced by the magnification bias, is much more 
pronounced (cf. Fig. 2b). For column densities beyond the peak location, N f(N) sharply 
cuts off. If the QSO is singly imaged, the ray-bending due to the lens increases the impact 
parameter of the QSO image compared to the case of inefficient lensing. If the QSO is 
multiply imaged, the brighter image whose spectrum dominates the absorption feature is 
also formed at a larger distance from the center of the lens than in the unlensed case. 
Thus, the lens prevents the light of the dominant image from passing through the highest 
column-density regions of the galactic disks, so that the highest column densities cannot 
be traversed by the dominant fraction of the QSO light. 

Since j/h > 0, equation (2.21) requires that 
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Figure 2.- Column-density distribution of damped Lya absorbers per logarithmic column-density interval, 
TV /(TV), for Ho = 50 kms -1 Mpc -1 . The integral of this function is proportional to the cosmological 
density of neutral hydrogen, J7hi- The heavy curve in both panels displays TV /(TV) without the influence 
of gravitational lensing. The three other curves were calculated for three different cosmological models, 
taking gravitational lensing into account. The cosmological model parameters are: J7o = 1, f2j\ = (solid 
line); i?o = 0.2, fl^ = (dotted line); and i?o = 0.2, fl^ = 0.8 (dashed line). In both panels, the QSO 
redshift is fixed at zqso = 3. In panel (a), the redshift interval for the absorbers is 2 < « a bs ^ 3, and in 
panel (b), 0.4 < « a bs ^ 1- While the influence of lensing is weak for high-redshift absorbers, a pronounced 
peak followed by a sharp cut-off arises for low-redshift absorbers. The QSO magnitude threshold is B = 16. 
We ignore the effect of possible obscuration by dust. 



Co < (Rh(N)) (3.1) 

for such lenses which contribute to damped Lya absorption with a column density > N. 
This implies that the maximum lensing efficiency of the absorbers is achieved for a column 
density iV pea k such that 

max Co = (i? H (iVpeak)) • (3.2) 

For N > iVp ea k, the redshift range for possible absorbers is clipped, so that their number 
decreases. The highest accessible column density iV cut is reached once 

min Co = (i? H (iVcut)) • (3.3) 

If min £o = 0, arbitrarily high column densities can appear, but only in an arbitrarily small 
number of absorbers. Therefore, the peak of f(N) at iV pea k is then followed by a decrease 
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with a finite slope for higher N . If min^o is finite, there exists a finite upper limit iV cut 
to the accessible column densities, where f(N) is steeply cut off. The minimum value of 
£o can vanish if and only if z\ =0 or Z2 = ^qso, as is the case in figure 2a. In all other 
cases, f(N) is cut off at a finite column density iV cut , as shown in figure 2b. Moreover, 
since £o oc cr 2 , the peak location measures the velocity dispersion and thus the depth of 
the gravitational potential of the absorber. 

In reality, the sharp peaks apparent in figure 2 may be broadened due to possible 
dumpiness of the HI distribution or due to an average over different density profiles of 
disk galaxies of different morphological types (see Broeils & van Woerden 1994). This 
smoothing of the peaks may render it difficult to constrain the value of the cosmological 
constant just from the measured location of the peaks. It is also important to note that for 
column densities ;> 10 22 cm -2 , the potential obscuration of the QSOs by dust is expected 
to be strong and therefore should be taken into account (Fall, Pei, & McMahon 1989; 
Pei, Fall, & Bechtold 1991; Fall & Pei 1993). The dust obscuration tends to reduce the 
observed number of QSOs with high column-density absorbers. Although this effect has 
the opposite sign to that of the magnification bias, its impact on the predicted f(N) is 
limited because of the sharp cut-off introduced by the bending of light rays due to lensing. 
We discuss this issue further in §4. We do not include a model for the dust obscuration 
in figure 2 because the central opacity of nearby disk galaxies is still a subject of current 
debate (e.g. Byun 1993; Davies et al. 1993; Byun, Freeman, & Kylafis 1994; Jansen et al. 
1994; Rix 1994), and the evolution of this opacity with redshift is even more uncertain. 

The distortion of the column-density distribution due to lensing affects the inferred 
value of the cosmological density parameter in neutral hydrogen, i?ni- Figure 3 shows the 
ratio between the inferred and the true value of f° r t wo choices for the absorber-redshift 
interval and for the three cosmological models specified before. Similar to observational 
studies (e.g., Lanzetta et al. 1995), we calculate the contribution to of damped 

absorption systems with 20 < log(iV cm 2 ) < 22. In figure 3a, we consider the absorber 
redshift interval < £ a b s < 1 while the QSO redshift ranges between 1 < £qso < 3. In 
figure 3b, the QSO redshift varies in the range < £qso < 3, and the absorber redshifts 
range within < z^ hs < zqso- 

The most prominent feature in both panels is the artificial decline of the inferred 
value of i?Hl with decreasing redshift, imitating a strong evolution of the neutral-hydrogen 
density with redshift. This noticeable effect arises because i?ni is dominated by the high 
column-density absorbers, for which gravitational lensing can strongly influence the result 
due to the magnification bias. While i?ni is generally overestimated for £qso i> 0.5, it 
can also be significantly underestimated for lower QSO redshifts. The reason for this 
underestimate is that the peak in f(N) shifts towards higher N for smaller QSO redshifts. 
This follows because the lensing length scale becomes smaller, requiring Rn(N) to be 
smaller and therefore N to be larger for high magnifications to arise. If the peak column 
density becomes larger than N2, the inferred i?ni rapidly drops below its true value (see 
Fig. 2). 

The discrepancies between the inferred and the true values of i?ni are strongest in 
case of a non-zero cosmological constant, as shown by the dashed curves. It is interesting 
that the artificial evolution in the apparent Qy{i{z) is already comparable to that observed 
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Figure 3.- The ratio between the inferred and the true cosmological density parameter in neutral hydro- 
gen, J7hi, is shown for two different choices of the absorber redshifts and the three cosmological models 
considered in figure 2. In the panel (a), we average the result over absorber redshifts in the interval 
^ ^abs ^ 1; while the QSO redshift is increased from zqso = 1 to zqso = 3. In panel (b), the QSO 
redshifts range within < zqso ^ 3, while the absorber redshifts lie between zero and zqso- In both 
cases, the inferred density parameter decreases at low redshifts. For low QSO redshifts, zqso < 0-5, the 
inferred density parameter underestimates the true HI density, but for zqso > 0-5, the neutral- hydrogen 
density is overestimated. A limiting QSO magnitude of B = 16 was assumed. For the sensitivity of the 
result to the QSO magnitude threshold, see figure 5. 



(White et al. 1993; Bahcall et al. 1993; Lanzetta et al. 1995), even before including any 
plausible amount of gas consumption by star formation (e.g., Kauffmann & Chariot 1994; 
Lanzetta et al. 1995). Although the observed evolution is often expressed in terms of i?ni 
as a function of z a b s rather than £qso, the previous searches were constrained to find high- 
redshift absorbers in the spectra of high-redshift QSOs (using ground-based observations), 
and bright low-redshift QSOs in surveys of low-redshift absorbers (using IUE data). 

The expected redshift distribution of the damped Lya absorption systems is presented 
in figure 4, given two choices of the absorber-redshift intervals. For the curves in figure 
4a, 2 < z a bs < 3, and for the curves in figure 4b, < £ a b s < 1- In both cases, the QSO 
redshift is kept fixed at £qso = 3. The solid curves show the result derived including 
lensing, while lensing was ignored for the dashed curves. Only results for the cosmological 
model with i?o =0.2 and Qa = are shown because for this model the difference between 
the lensed and unlensed cases is intermediate (cf. Fig. 3). Both panels suggest that it 
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Figure 4.- Normalized redshift distribution of 
damped Lya absorbers in the spectra of QSOs 
at a redshift zqso = 3, for absorbers at ei- 
ther high redshifts, 2 < « a bs ^ 3 (panel a) or 
low-redshifts, < « a bs ^ 1 (panel b). In both 
panels, the dashed curve was derived ignoring 
gravitational lensing, while the solid curves 
were derived including lensing. Only results 
for the cosmological parameters J7o = 0.2, 
Q A = are shown. The curves are arbitrarily 
normalized to unit area. 

might be difficult to distinguish between the lensed and unlensed cases on the basis of 
the absorber-redshift distribution alone. Especially if Qa is small, the curves are virtually 
indistinguishable given any reasonable error bars in the actual data. In addition, allowing 
for a range of QSO redshifts would tend to smear out and further reduce the difference 
between the two curves. 

The insensitivity of the absorber-redshift distribution to lensing (cf. Fig. 4), and the 
weakness of the distortion of f(N) due to lensing for high-redshift absorbers (cf. Fig. 2a), 
offer a plausible explanation as to why these imprints of lensing were not identified ob- 
servationally in the past. The lensing signatures would be most significant in samples of 
bright, high-redshift QSOs which are surveyed for low-redshift, damped Lya absorption 
with z a bs <J 1- Unfortunately, existing surveys of damped Lya systems (e.g., Lanzetta et 
al. 1995) do not include a sufficiently large number of cases with £qso J> 2-3 and z a b s <J 1. 
However, it should be straightforward to design a search strategy for the purpose of de- 
tecting the lensing signature, because high column-density absorbers and bright QSOs are 
both easy to detect. 





Figure 5.- The inferred neutral- hydrogen den- 
sity parameter depends on the limiting mag- 
nitude B of the QSO sample due to lensing. 



We consider a QSO redshift z 



QSO 



3 and ab- 



sorbers in the redshift interval < « a bs ^ 1- 
The three curves show the ratio between the 
inferred and the true fim for the three cos- 
mological models used in figure 2. 



Since the amplitude of the magnification bias depends on the QSO brightness, the 
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inferred value of depends on the magnitude limit of the QSO sample. The three curves 
in figure 5 show this dependence for the three cosmological models considered before. We 
assume £qso = 3 and < £ a b s < 1. The figure shows that the value of can be 
substantially overestimated for a bright QSO sample, but can also be underestimated if 
the QSO sample is faint. The latter effect occurs because the area in the source plane 
within which QSOs show damped Lya absorption with column density > N is reduced 
due to the focusing effect of the lens, as described before. The dependence of on the 
magnitude limit of the QSO sample can be used to separate evolutionary effects in the 
absorber population from effects caused by gravitational lensing. 
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Figure 6.- The fraction of multiple images 
among bright QSOs (B < 16) showing 
damped Lya absorption with N > 10 21 cm -2 
as a function of the maximum magnification 
ratio Hi / H2 ^ g of their images. In panel 
(a) we consider absorbers with 2 < « a bs ^ 3, 
while in panel (b), 0.4 < « a bs ^ 1- In both 
panels, «qso = 3. The different line types 
describe the three cosmological models con- 
sidered in figure 2. 



Figure 6 shows the fraction of all QSOs which show damped Lya absorption with N > 
10 21 cm -2 that are multiply imaged with a magnification ratio of their images (jii j [12) < <?• 
The two panels display results for different absorber-redshift intervals. For 2 < z a b s < 3, 
the fraction of high column-density absorbers that show double images can be used as a 
sensitive test for the cosmological constant. For low-redshift absorbers (0.4 < z a b s < 1), 
about 60%-80% of all QSOs should be doubly imaged with a magnification ratio below 
q = 2. The characteristic separation between the images is ~ 0.3" for an L* spiral galaxy 
at £ a bs ~ 1- Therefore, absorbers with N > 10 21 cm -2 can be used to flag possibly 
lensed quasars with small image separation. Since such absorbers are easy to identify 
through low-resolution spectroscopy, this approach can potentially enhance the efficiency 
of searches for new gravitational lenses with HST from its present value of ~ 1% (Bahcall 
et al. 1992; Maoz et al. 1993) up to several tens of percent. 



4. Conclusions 

The signature of gravitational lensing is particularly strong for systems that are easy to 
observe, namely high column-density absorbers (N ;> 10 21 cm -2 ) in the spectra of bright 
QSOs (B <J 18). There are five possibilities to verify the existence of lensing: 
1. The apparent cosmological density of hydrogen should decline considerably towards 
low redshifts (cf. Fig. 3). Such a decline has, in fact, been observed (White et al. 1993; 
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Bahcall et al. 1993; Lanzetta et al. 1995). However, the interpretation of the observed 
trend in Qy{i{z) is ambiguous because it could also result from the consumption of gas 
by star formation (e.g., Kauffmann & Chariot 1994; Lanzetta et al. 1995). To avoid 
confusion with such evolutionary effects of the damped-Lya absorber population, it is 
advantageous to make use of the dependences of the magnification bias on the source 
properties. 

2. The inferred value of i?ni should increase with increasing source luminosity (cf. Fig. 5) 
or source redshift (cf. Fig. 3) even if the absorber redshift interval is kept fixed. These 
effects cannot be imitated by an evolution of the population of damped-Lya absorbers. 

3. The column-density distribution for low-redshift absorbers in the spectra of high- 
redshift QSOs, /(iV), should show a pronounced peak followed by a sharp cut-off at 
high column-densities. For QSOs brighter than B ~ 17 with redshifts £qso J> 2, 
and for non-evolving absorbers in the redshift interval z a b s ~ 0.4-1, the peak should 
appear around N ~ (3-5) X 10 21 cm -2 . The location and height of the peak can 
be used to constrain the depth of the absorber potential wells or the value of the 
cosmological constant (cf. Fig. 2). Since i?ni is obtained through an integral over the 
column density distribution of damped systems, the first two signatures listed above 
may be easier to probe than this third one. This could explain why the first signature 
may have already been observed. 

4. When image splitting occurs, the absorption spectrum should show the superposition 
of two distinct absorption troughs corresponding to the different impact parameters 
of the two images. The two troughs are weighted by their corresponding image mag- 
nification factors (cf. Fig. 6). Since the equivalent width of each absorption feature 
scales as iV 1 / 2 , and since one of the images is usually faint, it may prove difficult to 
probe this spectroscopic signature. 

5. High-resolution imaging of absorbers with deep potential wells may uncover the exis- 
tence of multiple images. The image separation of ~ 0.3" (u c /220 kms -1 ) 2 expected 
for z a bs ~ 1 is easily detectable with HST unless the magnification ratio of the images 
is extreme. The fraction of doubly imaged QSOs that show damped Lya absorption 
with N > 10 21 cm" 2 ranges between 10% for 2 < z abs < 3 and 80% for 0.4 < z abs < 1. 
Therefore, strong damped Lya systems may serve as reliable flags for lensed QSOs 
with small image splitting (cf. Fig. 6). The identification of such systems requires 
only low-resolution spectroscopy. The selection of bright QSOs which are known to 
have strong damped Lya absorption lines may then be used to enhance the efficiency 
of searches for new gravitational lenses with HST (cf. Maoz et al. 1993). 

As in any other astrophysical study, there are various factors that may complicate the 
analysis of observational data. First, obscuration by dust may hide the lensed images of a 
QSO. The abundance of dust in damped Lya systems can be inferred from the reddening 
of the background QSOs. A comparison between the spectra of QSOs with and without 
damped Lya systems implies that the most likely dust-to-gas ratio in damped Lya systems 
at z ahs £ 2 is only 5-20% of the Milky Way value (Fall et al. 1989; Pei et al. 1991; Fall 
& Pei 1993). Therefore, the high-redshift absorbers are expected to be optically thin for 
HI column densities <J 10 22 cm -2 . Although the dust content of low-redshift absorbers 
could be higher than that, their obscuration is related to the low-extinction part of the 
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dust extinction curve in the V-band, where most QSO surveys are done. At low redshifts, 
the effect of dust can be readily corrected for by using the inferred obscuration profiles of 
nearby disk galaxies (e.g., Byun 1993; Byun et al. 1994; Jansen et al. 1994; Rix 1994). 

Note that the bending of light rays from the source due to lensing tends to weaken 
the effect of dust as it provides a minimum impact parameter, of order ~ 1-2 kpc, for the 
bright image of the QSO. Thus, if dust obscuration is strong in the inner part of damped 
Lya absorbers, then the significance of lensing may be emphasized even more, because 
without lensing fewer QSOs would have been seen behind the center of the absorber. 
Aside from adding QSOs through the magnification bias, the ray-bending due to lensing 
brings additional QSOs into view that would otherwise be hidden behind a thicker layer 
of dust. 

A second complication to our analysis involves possible dumpiness of the HI distri- 
bution in the absorbers, dumpiness would tend to broaden the distribution of column 
densities at a given impact parameter beyond the geometric spread induced by different 
inclinations. This effect would then broaden the peaks shown in figure 2. An argument 
against the existence of strong (;> 30%) dumpiness on small scales (<^ 1-2 kpc) is provided 
by the similarity between the column densities observed in the two images of the lensed 
QSO 0957 + 561 at £qso = 1.41; the separate lines of sight corresponding to the two images 
cross a common damped Lya absorber at z a b s = 1.391 and show a similar column density 
in their spectra (Turnshek & Bohlin 1993). In general, the higher the overdensity of the 
HI clumps, the smaller is their cross-section and the less likely they are to appear in ab- 
sorption spectra, dumpiness should be less important for estimates of averaged quantities 
such as than for differential quantities like f(N). 

Finally, a third complication to our modeling may arise from possible evolution in 
the properties of the absorbers. If the observed HI disks are much bigger at high redshift, 
then for a given HI column density the characteristic impact parameter increases and the 
lensing signatures are weakened relative to our no-evolution model. However, it should be 
easy to disentangle the pure lensing effect from any such evolutionary effect by examining 
the second signature mentioned above. 

Recently, there has been considerable interest in the cosmological abundance of 
damped Lya systems as a means of constraining variants of cold dark matter cosmolo- 
gies, such as the mixed dark-matter (MDM) model, in which galaxies form relatively late 
(Kauffmann & Chariot 1994; Ma & Bertschinger 1994; Mo & Miralda-Escude 1994). The 
magnification bias derived in this paper should in general be incorporated into those stud- 
ies. In addition, the observed sample of absorbers may be biased towards higher velocity 
dispersion a v because the magnification bias depends strongly on a v . However, this effect 
is limited by the exponential cut-off of the Schechter luminosity function. 

Some of the gravitational-lensing signatures listed above were not discovered observa- 
tionally in the past due to an unfortunate combination of selection effects. Ground-based 
observations were limited to absorbers at high redshift, while IUE observations (Lanzetta 
et al. 1995) focused on low redshift sources. In both cases the signature of lensing was 
artificially minimized by the lensing geometry. Forthcoming observations with HST could 
easily detect a few of the above five signatures if properly designed. The search strategy 
should focus on finding high column-density (N ;> 10 21 cm -2 ) absorbers in the spectra of 
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bright (B <J 16-17), high-redshift (^qso J> 2-3) QSOs. Detection of the peaks shown in 
figure 2 or the fraction of multiply imaged QSOs shown in figure 6 may be possible in a 
sample of more than ten damped Lya systems with N ;> 10 21 cm -2 and 0.5 <J z^bs <J 1- 
Quantitative measurement of these and the other lensing signatures can be used to con- 
strain the gravitational-potential depth of the absorbers and the cosmological constant. 
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